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Abstract 

The quantum discrete 4> 4 model at finite temperature is studied in 
the mean-field approximation. The phase diagrams are obtained for a 
wide range of the model parameters. The domains of applicability for 
the classical, quantum, and order-disorder limits are determined. It is 
shown that a remarkable part of the phase diagram is reproduced by the 
first quantum correction to the classical result. Parameters of Sn2P2S6, 
SrTiC>3, and KH2PO4 are placed in the phase diagram of the model. 

1 Introduction 

Quantum zero-point fluctuations play an important role in the ferroelectric phe- 
nomena for a wide range of ferroelectric materials (!]]. Some of the materials, 
such as SrTiC>3 , show a displacive (soft-mode) phase transition, while others are 
of the order-disorder type, for example the hydrogen-bond KH2PO4. Therefore, 
from the fundamental point of view, one should answer the two questions to 
characterize each material: (i) what is the interplay of quantum (zero-point) 
and classical (thermal) fluctuations and (ii) how close the transition is to the 
order-disorder or displacive scenarios. Working with model Hamiltonians, it is 
suitable to consider the quantum discrete 4 model here. The variation of the 
model parameters allows to go continuously from the displacive to the order- 
disorder second-order transition. An introduction of a finite mass and Plank 
constant gives rise to the zero-point fluctuations in the model. 

Although the discrete </> 4 lattice is one of the simplest possible microscopic 
models for the phase-transition phenomena, the studies have been restricted 
to the asymptotical cases so far. In the general case the transition type is 
neither purely displacive nor order-disorder and (or) fluctuations arc neither 
purely quantum nor thermal. Up to our knowledge, there is no a quantitative 
investigation of this situation so far. Moreover the exact values of the model 
parameters are not clear, than the known asymptotical cases realize. This paper 
is aimed, in particular, to solve the latter problem. 



Let us write down the explicit formulae and discuss the known asymptotical 
behaviour. The discrete c/> 4 model is a cubic lattice of 2-4 anharmonic oscillators, 
with a nearest-neighbour harmonic coupling. The Hamiltonian can be written 
as follows |]: 

« = + (i) 
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Here x n is the displacement of the n-th atom in the unit cell, A, B, C are pa- 
rameters of the model. The first two terms in the potential energy V form a 
double-well potential and the last term stands for the interaction. There is only 
the inteaction between nearest- neighbours in the model, i.e. a nn i — 1 if n,n' 
are the nearest neighbours and vanishes elsewhere. Note that the Hamiltonian 
(1) includes the displacement of atoms, but not the actual distance between 
them. 

If there were no any fluctuations in the system, all atoms would occupy 
the same position in the cell x n = sjAjB (or x n = —y/A/B). There are two 
sources of fluctuations. The dimensionless mass m — MCA 2 /h 2 B 2 determines 
the role of quantum zero-point fluctuations in the model (the system is classical 
at m — > oo). The classical thermal fluctuations are governed by the reduced 
temperature t = TB/AC (the system becomes purely quantum at t = 0). Large 
fluctuations can destroy the ordering, so that the average of x vanishes. The 
system of a dimensionality d > 1 shows a second-order transition for any pos- 
itive A, B,C. The parameter a = A/C determines the type of the transition: 
displacive transition corresponds to a — > 0, while the limit a — > +oo corresponds 
to the order-disorder system || Qj. At a given a there is a critical line t c (m c ) 
in the t — m phase diagram (hereafter the subscript "c" stands for the critical 
values of temperature and mass). 

In the reduced dimensionless variables, the Hamiltonian takes the form 
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The asymptotical cases are studied extensively. Basically, they are the fol- 
lows. 

At the limit a — > +0 the transition is of the soft-mode type. This asymptotic 
is studied at t = 0, as well as at finite temperature || ||. Since the double- well 
potential is small for this case, it can be treated as a perturbation. Therefore 
a convenient method of the investigation of this limit is the independent-mode 
approach ||. The system is harmonic in the zero approximation, therefore the 
phonon vibrations Yfe oc e Xj can be treated as uncoupled oscillators so that 
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one assumes < YkYi >oc 8(k + I) (hereafter the triangle brackets stand for the 
averaging over thermal and quantum fluctuations). In the first approximation 
the anharmonic term renormalize the phonon dispersion. Virtually, the real 
interaction can be reduced to an effective (i.e. dependent on the temperature 
and mass) harmonic one. At certain values of t and m the zero-frequency 
mode (soft-mode) appears in the system indicating the phase transition. The 
independent-mode approach is exact in the a — > +0 limit. However, serious 
difficulties appear for a finite a; particularly, t c and m c are incorrectly predicted 
to be independent on a in the independent-mode scheme. 

At the opposite limit a — > +00, the discrete </> 4 model shows an order-disorder 
phase transition. The physical picture in this case is the follows. Since an on- 
site potential is much larger than all other energy scales, atoms are localized 
in either left or in the right side of the double-well in the zero approximation: 
x n w ±yj A/B. Therefore, each oscillator can be treated as a two- level system. 
The tunnel coupling between the two sides of a well produces some energy 
splitting r, which depends on the mass m implicitly. This case coincides with 
the transverse-field Ising model [M @j. This model describes an array of \- 
spin particles in the external field with a nearest-neighbour interaction between 
particles. The Hamiltonian of this model is 

i " i,j 

where Sf, Sf are | spin operators. The phase transition can be observed 
by changing the two parameters: temperature t and transverse field T. The 
transverse-field Ising model is also widely studied at t — and at finite temper- 
ature by various approximations. High-temperature series expansions, ground- 
state perturbation theory, correlated-basis-function analysis, random phase ap- 
proximation and its generalizations, 1/3? and 1/5 expansions, continuum field 
theory methods have been used @, |, |, [l0[ ||[ |3| . 

The classical case corresponds to the condition m — > 00. In this case the 
momentum-coordinate uncertainty is small, so that the partition function can be 
factorised as Z = J exp(— ^p^/2mt)d Ar p„ ■ J exp(— v/t)d N x. The temperature 

n 

dependence of the order parameter is determined here by the second term and 
is therefore independent of mass m. The dependence t c (a) is studied both 
numerically and by analytical approaches for 2D, layered and 3D systems ||. 
An extension of the mean-field scheme has been developed, which is appropriate 
for the description of Monte Carlo results with a good accuracy in a wide rage 
of parameter a. 

The quantum limit takes place at t = and at finite mass m. The de- 
pendence m c (a) is studied jl4| by quantum Monte Carlo simulations and the 
simplest analytics for a wide range of a. 

In this paper, we study the phase diagram of the quantum discrete <fi model 
for a full range of parameters, and determine the domain of applicability for 
the asymptotical results mentioned above. We also consider the leading-order 
quantum correction to the classical result. In a conclusive part of the paper the 
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parameters of several ferroelectric materials are placed in the phase diagram 
obtained. 

We chose the mean-field (MF) approximation for the study because of its 
simplicity. Actually, we are interested in quite crude characteristics such as the 
domains of applicability. Therefore we believe that this method is suitable here, 
although it is not very accurate. The validity of the MF approach is discussed 
in more details in the end of the paper. 



2 The quantum mean-field approximation 

In the MF approach, the real interaction between particles is replaced with an 
interaction of the on-site oscillator with an external field E — Ad < x >. The 
self-consistent set of MF equations for < x > is formed by this formula and the 
expression for the average displacement of the on-site oscillator [jl5| 

_ Sp ( xe -{H Q -E<x»/T^ 

<X> ~ Sp(e-( H »- E <*»/ T ) ' (4) 

where 

Eo = ^ + V , (5) 

V = LlC-fjat + jX*. 

In this paper, we are interested in the position of the transition point in phase 
diagram. The condition < x >— > is fulfilled in the point of the phase transi- 
tion. Keeping the leading-order in < x > terms, one obtains the MF equations 
for the critical point: 

AdCxi = 1, (6) 

where xi 1S the static linear susceptibility, which is determined at the finite 
temperature by the standard quantum-mechanical formula Jl6[ : 

E2dy F — Ei 

tj 

Here, F = —TlnZ is the free energy (Z is a partition function), Ei =< i\Ho\i > 
and dij —< i\d\j > are respectively the eigenenergies and the dipole matrix 
elements for the on-site Hamiltonian Ho, defined in Eq. (j^). 

Let us briefly discuss the structure of the spectrum. In the displacive limit 
a — > +0 it is almost equidistant, because the nonlinearity of the potential is 
small. In the order-disorder case of a — > +oo the situation is quite different. 
There are several doublets of the lowest levels, whereas the rest of the spec- 
trum is situated much higher (Fig. 1). A number of doublets increases with 
an increase of m, since the system becomes more classical and the spectrum 
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becomes almost continuous. At the same time, the behaviour of the system is 
determined by the first doublet at low and moderate masses. This is the case 
of the transverse-field Ising model (3). 

The spectrum and matrix elements of the on-site oscillator (^) were cal- 
culated numerically. Here are some details of the numerical procedure. The 
eigenfunctions of the Schrodinger equation with Hq were found by the phase 
method p7| . We controlled the accuracy of calculations by the check of the 
condition < i\j >= Sij, which was fulfilled with an error bar not exceeding 
10~ 3 . We used approximately 50-70 levels of H to obtain the phase diagrams 
m c (t c ) in the quantum MF approach at moderate values of a and for m « 5. 
This amount is necessary to describe correctly the case of large mass, than the 
levels are quite dense. The number of necessary levels is increased with an 
increase of t as well. 

The phase diagrams are obtained for a wide range of the parameter a. They 
correspond to the displacive limit, order-disorder limit and to an intermediate 
case. The dependencies m c (t c ) are plotted in Fig. 2 by solid lines for the three 
values of the parameter a for 2D and 3D lattices. All curves show a good 
agreement with the known asymptotes for the quantum limit (t = 0) and for 
the classical case (m — » +oo) fi ||, Q. 

3 The limits 

Now we turn to the discussion of the domains of applicability for the four men- 
tioned asymptotical cases: displacive, order-disorder, classical and quantum 
limits. In the displacive limit, the MF scheme suffers serious difficulties (see 
discussion in the end of the paper). However, from the previous studies of the 
classical g, || and quantum Q systems, we can estimate that the displacive 
limit occurs at a smaller than awl. For the three other limits, we use the 
following scheme. We calculate the critical parameters by the MF expressions 
which are valid in the corresponding asymptotical cases, and compare them with 
the direct MF calculations. For each limit we determine the areas where the 
error bar does not exceed 1% and 5%. 

For the quantum limit a deviation of m c from the value of m c (t = 0) is 
detected for each a. The two areas, where m c differs from m c (t = 0) in a value 
smaller than 1% and 5%, respectively, are marked in Fig. 2. 

To determine the areas, where the system can be considered as classical, 
we compare t c (m c ) with the classical value i c (oo), at a given value of a. The 
obtained areas, where the classical result delivers a 1% and 5% accuracy, are 
shown in Fig. 3. For comparison, we also present the "purely quantum" MF 
dependence of m c {a) at t = in the same figure. 

Finally, the domain of applicability of the order-disorder limit is estimated 
from the comparison of the dependencies m c {t c ), obtained from the calculation 
by (0) with a complete set of eigenfunctions, and with only the first doublet in 
the spectrum. The result is shown in Fig. 3. 

Thus, the values of the parameter a, which are used to plot the curves in 
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Fig. 2, cover a complete range of phase transitions from the displacive type for 
small values of a, to the order-disorder one for large values of a. Besides, the 
presented intervals of m and t (Figs. 2, 3) show a crossover between quantum and 
classical phase transitions for the discrete </> 4 model. 



4 The quantum correction formula 

We have also plotted the leading-order quantum correction to the classical MF 
value. The idea is to use usual classical MF scheme, which takes into account 
quantum fluctuations in the system. It can be obtained in Fcynman's con- 
sideration of the path- integrate representation for the partition function ]l9| . 
The partition function of the on-site oscillator Z = Sp (er( H °~ E<x> >/ T \ in this 
representation can be written as follows: 

Z = J [Dx] exp \h~ l T L{x{T))d^j , (8) 

where the integration is taken over all classical trajectories x(r), describing 
the propagation in the imaginary time r = 0...H/T with periodic boundary 
conditions applied. L(x(t)) is a Lagrangian along the trajectory x(t). The 
semi-classical approximation consists in the assumption that the deviations of 
the trajectory from its mean position xq = T%~ x J xdr are small. One should 
perform a Taylor-series expansion of (||) in powers of x(t) = x(t) — xq and keep 
the terms up to x 2 . After this, x can be integrated out, yielding 



Z « Z Q J dx exp -i (v (x ) - Ex 



h 2 8 2 Vr 



2AMT dx 2 J 



(9) 



Here Z is a factor, independent of Vq(x). Thus, in this semi-classical approx- 
imation we can consider a quantum system using its classical potential with 
quantum correction pSj : 

Vl -V + —^. (10) 
q ° 24MT V 1 

Virtually quantum correction renormalises the on-site potential parameters 
making them dependent on temperature. 

The results obtained from the classical MF theory with the quantum correc- 
tion ( |Io| ) are presented in Fig. 2 by the dashed lines. All curves demonstrate an 
excellent agreement with ones obtained directly from the quantum MF calcula- 
tions (^). Certainly, the approach ( |l0| ) does not work at low t and m at all. In 
our case of discrete <j> model we can calculate the phase diagrams in this way 
down to the mass m ~ 0.35. Nevertheless, the area where the semi-classical 
approximation works is large enough to demonstrate a crossover between the 
quantum and classical phase transitions. We would like also to stress that the 
calculation via the quantum-mechanical formula (|^) is about 10 4 — 10 5 times 
longer than the classical MF method with the quantum correction (|l0|). 
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5 Experimental data and discussion 



The estimation of the parameters m and a of the discrete </> 4 model for real 
materials and their accordance to the limit cases on the phase diagram m c (a) 
(Fig. 3) is of interest. We will discuss a situation with the following ferroelectrics: 
Sn 2 P 2 S 6 , SrTi0 3 , KH 2 P0 4 (KDP). 

The transition in Sn 2 P 2 Sg is of an intermediate type |2(], |21j| . Therefore, 
there are thermal regions in phase diagrams, where Landau theory is valid, i.e. 
the free energy is of the form: y « a ^ P — h P( T J P , with a oc (T — To), 
/3 « const. In our estimations for Sn 2 P 2 Sg, we use the data for the Curie 
low and the value of the saturated spontaneous polarization to obtain values of 
coefficients a(0) and /3(0) (^2). The following estimation for the parameters m 
and a can be suggested (for an estimation of a, our method is very similar to 
what was used in ^2|). The following relation can be easily obtained: A 2 /B = 
V uc ct 2 1 (3 (V uc is a volume of the unit cell). Also, accordingly to the phase 
diagram for the classical discrete tfi 4 model 0], one has kT c = XAC/B, with 
A « 6 for an intermediate type of the transition. From these two relations one 
can guess about the ratio a = A/C « \V uc a 2 / kT c (3 . A suitable estimation for 
the parameter m = MCA 2 /h 2 B 2 comes from the observation that the mean- 
square of the displacement < Air 2 >« A/B is known from the experiment pTj| , 
p3[ . Combining this with the two previous formulae, we obtain m ss MkT c < 
Ax 2 > /h 2 \. The values a(0), /3(0), V uc , kT c and < Ax 2 > obtained from the 
experimental data [ p0| , |2l| , p2| , p3| , as well as our estimations for m and a are 
presented in Table 1. 

We also use the same method in estimations for SrTi03, which is a displacive- 
type material However, there are some special issues in this case. An 
ordering in SrTiOa is absent due to the quantum effects. The dependencies a(T) 
and P(T) are Landau-like only at temperature down to approximately 50 K; at 
lower temperature they are saturated. Therefore, in our estimations we use the 
experimental data ||, || for a{T),(3{T) for T = 50..100K. We approximate 
these data by the Landau-like dependencies a oc (T — T c ), f3 — const, i.e. 
assuming a "classical" picture of the phase transition. We extrapolate these 
"classical" dependencies down to T = 0, and obtain the "classical" critical 
temperature, and a(0),/?(0) from this extrapolation. If the system were purely 
classical, a(0),/3(0) and T c would take these values. The extracted data are 
presented in Table 1. Since the displacive phase transition takes place in SrTiO,3 
we use the estimation for A as A ~ 3, according to the previous study of the 3D 
discrete cffi model All other values required for estimations of a and m are 
known from the literature [^4| ^6) and arc placed in Table 1 as well. 

The corresponding points for Sn 2 P 2 Sg and SrTi03 lie in the phase diagram 
(Fig. 3) in the classical and quantum areas, with the intermediate and displacive 
types of the transitions, respectively. These data correspond to the experimental 
observations |,00gJ|||3[||§6|]. 

However, the Landau theory is hardly applicable even qualitatively in the 
case of KDP due to the order-disorder phase transition of the essentially quan- 
tum nature (the transverse- field Ising model is quite suitable here). We make 
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Table 1: Experimental data and estimations for m and a for Sn2P2Sg and 
SrTi0 3 . 

a(0), J-m-CT 2 /3(0), J-m 5 -C- 4 V uc , m 3 kT c , J < Ax 2 >, m 2 m a 

Sn 2 P 2 S 6 2.7 x 10 9 10 11 2.3 x 1CT 28 4.6 x 10~ 21 10~ 21 12 21 

SrTi0 3 5 x 10 7 4 x 10 9 6.4 x 10~ 29 7 x 1(T 22 1.7 x 10~ 22 0.3 0.2 



use of the following method of the estimation of parameters. The splitting of 
the first doublet of the on-site oscillator spectrum |27) is known to be about 
lji — luq 80 K, while the third level lies on the distance of approximately 
u)2 — ujq ~ 800 K from the first one. Further, we take into account the value for 
T c w 120 K. Given these data, we adjust a and m c to deliver the known values 
of (toi — luq) I (102 — wo) and [lo\ — luq) /T c . This estimation gives us a w 250 and 
m ~ 0.1. Thus, the point for KDP disposes in the quantum area of the phase 
diagram (Fig. 3) and corresponds to an order-disorder transition. 

Thus, our estimations place the three mentioned materials to the asymp- 
totical cases they are commonly attributed to. Now we discuss how the points 
can move on the phase diagram. The most direct way is to change the isotope 
mass of the proper components. An elegant example is the restoration of the 
ordering via the change from 16 O to 18 O in SrTiC>3 p8fl . The transition tem- 
perature happens to be about 23 K in this case, while the "classical estimation" 
for T c (placed in Table 1) is approximately two times larger. This corresponds 
well with our data: at small a a 5-7 % change in m can increase the transition 
temperature from to t c « 2, while the "classical" asymptote is t c m 4 (see 
Fig. 2). Clearly, however, the "isotopic" change in m cannot be large. 

Finally, we discuss the applicability of the MF method for the present prob- 
lem. As the MF scheme ignores correlation effects in the system, it overestimates 
the transition temperature (or underestimates the critical mass). This overesti- 
mation is particularly large for the displacive limit. Namely, the MF approach 
predicts a finite value of t c (a — > +0) in any dimension, whereas actually fluctu- 
ations result in t c cx In -1 a in 2D at finite temperature ^, [l8|]. However, for the 
moderate and large values of a in 2D, as well as for the whole range of a in 3D, 
we expect that the qualitative behaviour of t c (a) and m c (a) is reproduced cor- 
rectly. Quantitatively, it is known that the MF scheme overestimates t c |j, or 
underestimates m c [[u] by some 30-50 %. We expect that the similar situation 
takes place for the general case, then both thermal and quantum fluctuations 
are present. To get more justification, we have now started the quantum Monte 
Carlo simulations for the system under consideration . The pilot calculations 
support the qualitative validity of the MF results. 
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Figure 1: The on-site potential and its spectrum for the two limits of the discrete 
(ft 4 model. The left panel is a typical spectrum for the order-disorder limit 
(a = 100). The right panel is a typical spectrum for the displacive limit (a = 1). 
Parameters m = 1 and d — 3 are used for the both cases. 




Figure 2: The dependencies of the critical mass m c on the critical temperature t c 
for the three values of the parameter a: a — 4, a = 16, a = 64. The logarithmic 
scale for m c is used. The left and the right panels show m c (t c ) for 2D and 
3D, respectively. The solid lines show the result of quantum MF approach (Q), 
the vertical lines are the classical MF limits || [I| the dashed lines are the first 
quantum correction (10). The lines marked as 1% and 5% show the estimations 
for the domains, where the "quantum limit" value of m c (t — 0) delivers the 
corresponding accuracy for m c (t). 
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Figure 3: The phase diagram for the limit domains of 2D and 3D systems (the 
left and the right panel, respectively). Solid lines indicate the classical limit 
m — > +00 (hatching on these curves indicates the beginning of the classical 
area with an accuracy of 1% and 5%). The purely quantum limit corresponds 
to the curve with t = (the domain of applicability for this limit is shown in 
Fig. 2). The dashed lines indicate the order-disorder limit area. Vertical lines are 
the classical asymptotes for these curves. The displacive limit approximately 
corresponds to a < 1. Points on the right panel are the estimations for m and 
a from the known experimental data for Sn2P2Sg, SrTi03, and KH2PO4. 
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